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Abstract 

The canonical ADM equations are solved in terms of the conformal factor in the 
instantaneous York gauge. A simple derivation is given for the solution of the two 
body problem. A geometrical characterization is given for the apparent singularities 
occurring in the AA-body problem and it is shown how the Gamier hamiltonian system 
arises in the ADM treatment by considering the time development of the conformal 
factor at the locations where the extrinsic curvature tensor vanishes. The equations of 
motion for the position of the particles and of the apparent singularities and also the 
time dependence of the linear residues at such singularities are given by the transfor- 
mation induced by an energy momentum tensor of a conformal Liouville theory. Such 
an equation encodes completely the dynamics of the system. 

1 Introduction 

The problem of providing solutions to 2 + 1 dimensional gravity [0 has been approached 
along different lines. In absence of particles the dynamics becomes non trivial only for closed 
universes of genus higher or equal to one. Moncrief || and Hosoya and Nakao || gave the 
hamiltonian treatment of this problem providing the complete reduction of the hamiltonian 
to the physical parameters i.e. the moduli of the time slices. This procedure was obtained 
in the York K = t gauge; in the case of the torus the problem can be dealt with explicitly 
while for higher genus, even if the problem is well defined we do not posses an explicit form 
for the reduced hamiltonian. 

A different approach was put forward by 't Hooft Q by describing the evolving Cauchy 
surfaces in terms of polygonal tiles which join along segments where the extrinsic curvature 
is singular. The dynamics of the system is codified in transition rules which intervene when 
e.g. the length of a side of a polygon goes to zero or when a particle collides with a side of 
a polygon. Such an approach is applicable both in presence and in absence of particle and 
it is suitable for an algorithmic or numeric approach. A similar but different approach was 
given by Waelbroek p. 
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In presence of particles progress was made in the papers by Bellini, Ciafaloni and Val- 
tancoli || and by Welling in the first order formalism by going over to the instantaneous 
(K = 0) gauge; in these papers it was shown that the problem is equivalent to a Riemann- 
Hilbert problem. It can be solved explicitly for the case of two particles, in terms of hy- 
pergeometric functions. For three or more particles one encounters a feature known in the 
mathematical literature as apparent singularities. 

In this paper we shall consider 2 + 1 dimensional gravity in presence of massive particles 
by exploiting the hamiltonian formulation. Thus our setting is a second order formalism and 
the basic equations are those derived systematically from the variation of the ADM action. 
The reason for such a development is to give a treatment which resides completely within 
the canonical framework. All results obtained in [J7| come out in simple fashion from such 
an approach; moreover we shall prove that the Gamier hamiltonian system for the apparent 
singularities is a direct outcome of the canonical ADM formalism. Again the exploitation of 
the instantaneous York gauge plays a major role; the technical advantage of such a gauge 
is to reduce an equation of sinh-Gordon to one of Liouville type to which powerful methods 
of complex analysis apply. Unfortunately the applicability of such a gauge is restricted to 
open universes. 

We shall see that the described approach provides an elementary way to solve the two 
body problem: the exact solution of the motion can be derived even before the explicit 
computation of the metric; combined with the qualitative knowledge of the asymptotic 
metric it gives a complete description of the scattering. The exact metric is obtained by 
solving a Liouville equation, leading of course to the same results obtained in || and @. 

One of the aims of presenting a canonical solution of the problem is to provide a framework 
for the quantization of the theory even if here we shall confine ourselves to the classical 
problem. 

Quantization schemes for 2 + 1 dimensional gravity in absence of particles have been 



proposed in H, EJ, [10[ and in presence of particles in [O, [L2|, O] . 



The key quantity of our approach will be the conformal factor which describes the metric 
of the time slices. The lapse function N is given by the derivative of such conformal factor 
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with respect to the total energy and the shift functions are given in terms of derivatives of 
the lapse function. 

As mentioned above, with more than two particles there are apparent singularities and 
one has to provide evolution equations for their positions and residues. We show that all this 
information in contained in the ADM equation for the time derivative of the conformal factor. 
In fact by exploiting Schwarz's relation between the coefficient of the fuchsian differential 
equation and the reduced conformal factor one obtains the Gamier equation by equating the 
residues at the polar singularities at the position of the apparent singularities. 

In the present paper we shall outline the main features of the approach; technical details 



are found in ref. 14 



2 The action 

The action for the gravitational field including the boundary terms is given by |15| 

~ SGrav = [ dt d 2 Xy/g 11+ 

c 4 Jm 
2 / El d 2 xJgK + 2 / d 2 xy/^B+ (1) 

JS JB 

2 / dxy/a r] 

where g^, TZ are the 2 + 1 dimensional metric and curvature, the 2-dimensional metric 
of the constant time slices E t ; S and Ei the initial and final time slices; the second 
fundamental form of the time slices E t and 0^ the second fundamental form of the lateral 
boundary B whose volume form is ^—7; is the future pointing unit normal to the time 
slices and the outward pointing unit normal to B; B t = S t R B, yfa is the volume form 
induced on B t and sinhr/ = n^. The boundary terms are going to play a fundamental role 
in the computation of the reduced action; this paper however will be based completely on 
the equations of motion which do not depend on the explicit form of the boundary terms. 
In presence of particles we must add to ([[]) the term 



m r 



-^u{q n )qnq u n (2) 
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with m n = 167r( gv m " ; o < m n < An. In the ADM notation where the metric is written as [IB 



ds 2 = -N 2 dt 2 + g^dx 1 + N i dt){dx j + N j dt) (3) 
the volume term apart for a divergence, can be rewritten in terms of the canonical variables 



as 



300 

/ dt d 2 x W j gij - WHi - NH] (4) 

J M L J 

where we have introduced the momenta canonically conjugate to g^. Similarly the matter 
action takes the form 

fdtY J {Pmq i n + N\qn)Pni- (5) 



N(q n )^/P ni P nj g^(q n )+m 2 n ). 

3 The constraints 

In (2 + 1)— dimensions, the great advantage of a York-type gauge, d{K = 0, is to ensure, when 
combined with the conformal gauge g^ = exp(2cr)5 i j, the decoupling of the traceless part 
7r Tu of the momentum from the spatial metric in the diffeomorphism constraints^. Indeed, 
the diffeomorphism constraints, 

Hi = -2^7,^ - - ^)Pm = 0, (6) 

V 9 n 



become 



and 



d z 7r Tz - z = -W5 2 (x-q n )P n 



d,n T " z = -W5 2 (x-q n )P n 



which are solved by 



7T 



1 P 

Tz _ nz 



i-Ezz^ + P'M (7) 



5 This decoupling reflects the fact that the divergence of a traceless symmetric tensorial density of rank 
two is Weyl-invariant in two dimensions. 
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and its complex conjugate expression. Here we introduced the complex coordinates z and 

z with z = x + iy and P z z (z) is an entire function. In the following, we shall adopt the 

maximally slicing gauge K = 0, where coincides with its traceless part so that so that 

eq.(^) holds for tt z z and its complex conjugate for n z z . 

We come now to the Hamiltonian constraint, 

1 



71 a - c) 



H= — 

+ Y / S\z- z n )^ml + P ni gvP nj = 0. 

n 

If we integrate the l.h.s of eq.([8D over a slice of constant time, an elementary application 
of the Gauss-Bonnet theorem shows that the condition K = is consistent only with open 
universes or with closed universes of genus 0. In the latter case consistency can be achieved 
only for static universes as it will be shown in Sect. 4. Moreover, for open universes, eq.(J|) 
determines the asymptotic behavior of the momenta 7r* • in terms of that of a. In fact if 
we assume that the asymptotic metric has to be conical, e 2cr ~ {zz)~^, /i > 0, the above 
equation implies that 

J d 2 z 27r 2 f 7r 5 2 exp(— 2er) < oo. 

This, in turn, entails that n z z at infinity vanishes faster than 1/z. Thus, the requirement 
of having a conical metric at infinity sets the function P z z {z) to zero and in addition it 
imposes that J2 n Pnz = 0. Because of the explicit form of the solution (^), we actually end 
up with a tt z z that goes to zero as 1/z 2 . As a consequence, if we define, mimicking the four 
dimensional counterpart, the total momentum P to t of the system as the surface integral of 
71"^-, we are led to conclude that P tot vanishes identically. It is worth noticing that in the 
present gauge P tot is simply equal to J2 n Pn an d no gravitational correction appears. The 
clash in (2 + 1)— dimensions between the boundary conditions ensuring finite energy and 
those ensuring finite momentum was first noticed in |Tj| |1^, |2(J. In particular, in [PJ the 
origin of this conflict was traced back to the impossibility of performing asymptotic boosts 
without introducing spurious singularities in the metric. 

We return now to eq.@. It is not difficult to show [14| that e~ 2cr vanishes on the 
particle singularities and as a result the hamiltonian constraint assumes the form, after 
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setting 2a = 2a - \n(2n~ z z ii z 2 ) 

2Aa = (9) 

B n 

The above equation is an inhomogeneous Liouville equation. The sources are given by the 
particles and by some other auxiliary sources whose position zb are the zeros of ir z z and their 
strength is fixed to —Air. The conformal factor a is actually regular at the points zb and for 
this reason we shall refer to them as apparent singularities. 

We remark that the reduction of the hamiltonian constraint to a Liouville equation is 
the result of the K = condition. The K = const 7^ gauge gives rise, as seen from eq.(^]) 
to an equation of the sinh-Gordon family which is of much more complex nature. 

Finally we note that the conformal gauge i.e. gij = e 2cr 5ij is still not a complete gauge 
fixing. One is allowed to perform analytic transformation on z which leave, in the case of 
open universe, the point at infinity fixed i.e. time dependent translations, rotations and 
dilatations. This freedom can be exploited e.g. to keep one particle fixed at z\ = and 
an other fixed at Z2 = 1. In general such a fixing is equivalent to a frame which rotates at 
infinity. 



4 Equations for the lapse and shift functions 

Before discussing the solutions of the inhomogeneous Liouville equation we shall write down 
the equations for iV and N l ] these are an outcome of requiring that the gauge conditions are 
preserved in the Hamiltonian evolution. Combining the canonical equation for 7T*- 7 and the 
one for the relation K = (or equivalently n l i = 0) takes the form 

AN = e- 2a N. (10) 



Again we have exploited the vanishing of e at the particle sources [14]. Since we have 



chosen the conformal gauge, the traceless part of the equation for determines N\ i.e. 

d- z N z = -n z - e~ 2a N (11) 
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and its complex conjugate. These complete the list of our equations for the metric. By 
applying Gauss theorem to eq. ([T0|) one easily derives the asymptotic behavior of iV for large 

\z\ 

N ~ ln(zz). (12) 

Such an asymptotic behavior of the lapse function is typical of the instantaneous York gauge 
and should be contrasted with the behavior N — > const, of the familiar DJH gauge IB]. The 
behavior (0) holds both when the plane describes the time slice of an open universe and 
when it describes the sphere; but in the latter case in order for the scalar N to be well 
defined on the Riemann sphere is should go to a constant at infinity. This occurs only for 
7r l j = i.e. when the equation for N reduces to AN = 0. Thus we see that except for the 
trivial stationary case in which P n = the York instantaneous gauge applies only to open 
universes. 



5 The Liouville equation 



As is well known, the problem of solving the inhomogeneous Liouville equation ([5]) can 
be reduced to that of solving a fuchsian second order differential equation of projectively 



canonical type EBJ 



2/'(C) + Q(CMC) = o 



where 



Q(z) = E 



+E 



B I 



1-/4 , Pn 

4(z - z n ) 2 2(z - z, 
3 (3 B 



(13) 
(14) 



+ 



A(z - z B f 2{z - z B ) 
If we define the function f(z) as the ratio of two independent solutions of eq. fll3|) , 

kyi(z) 



(15) 
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the solution of eq. (|9]) is built from the Poincare metric for the pseudosphere by performing 
on it the conformal diffeomorphism defined by eq. fll5D 

a - mmr ut> ' 

The poles in ( |I4"D include both particle and apparent singularities. The residues of the second 
order poles in Q(z) are directly related to the strength of the delta singularities in eq.©: 
in fact, /i n = m n /47r while the second order residue at the apparent singularities is fixed 
to be —3/4. The latter value is required by the regularity of 2a at Zb- The determination 
of the first order residues is more difficult and it is related to the knowledge of the explicit 
solution of the Riemann-Hilbert problem associated to flTBl). They are fixed in principle 
by the requirement that the conformal factor a be monodromic in the complex plane. As 
eq.([l|]) can be rewritten in the form 

e - = - 7^" (17) 
(2/22/2 - kk VxyiY 

where w\2 is the constant wronskian, it is seen that such a requirement is equivalent to the 
imposition that the vector (ky 1 ,y 2 ) transforms under a representation of £77(1, 1) when one 
encircles a singularity in the complex plane. It provides the relation between the geometric 
structure approach |22| and the one here described. 

An important feature in the solution of eq.(|9]) is that one has to give the residue of the 
singularity at infinity and this information is not contained in the parameters appearing 
in the sources. In different words, in solving eq.([|) the total energy of the system has to 
be specified as an arbitrary parameter M because the positions and the momenta of the 
particles at the initial time are not enough to determine its value. The above conclusion 
might seem counterintuitive especially if we compare it with what happens in the DJH 
geometrical approach, where the knowledge of all z t and Pi fixes M completely through a 
nonlinear composition rule for the momenta. This apparent mismatch simply reflects their 
choice to have, in contrast with the present approach P, |7|, IT3}], a flat (not globally defined) 
metric. As a consequence, in their approach, all the dynamical information is shifted into 
the particle variables. In our approach a complete set of initial conditions for the conformal 
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factor a is instead provided by any unconstrained parameterization^ of the rational function 
Q(z): an example will be given in sec. 9. 

The conformal factor 2a plays the key role in all subsequent developments; once the 
solution of eq. (|jD is known the equations for N and N l are solved as follows. As the sources 
in eq.(|9|) do not depend on the total energy M, by taking the derivative with respect to 
M we find a solution to the homogeneous equation fllPf ) for N i.e. apart for an arbitrary 
normalization constant we have 

As for the equation ( |TT|) multiplying it by 2n z z and using eq. (|T0|) and taking into account 
that n z z is a function of z we get 

N z = -^--d z N + g(z) (19) 

where g(z) is a function of z. A similar solution holds for N z . The function g(z) has to be 
chosen as to cancel the poles generated by the zeros of n z z (z), which occur only in presence 
of three or more particles, and if we are interested in describing a reference system which 
does not rotate at infinity g(z) has to be chosen as to give N z < \z\ at infinity. We shall see 
in sect. 7 that, due to J2 n Pnz = 0, such a function can diverge at infinity at most linearly. 
As we shall see g(z) encodes the time evolution of all the quantities of the problem. 



6 The particle equations of motion 

The canonical equations for the particle positions z n and momenta P n are obtained by varying 
the action with respect to z n and P n and are given by 

z n = -N i (z n ,t) = -g(z n ), (20) 

6 In fact the j3 n and /3 a are not all independent: In the next section we shall see that fuchsianity and 
the value of the residue at infinity impose two sum rule on them. Moreover, the requirement that conformal 
factor a be regular at za fixes the Pa in terms of the other parameters by means of the so-called no-log 
condition, see sec. 9. 
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nz na ~dz~ ^ n ~dz~ ) 
and their complex conjugate expressions. We note that z n is simply given in terms of the 
meromorphic function g(z) as in eq.([H]) ir z z diverges at the particle positions while d z N stays 
finite. It is of interest that the r.h.s of eq.([H|) depends only on the function g(z) and on the 
linear residue (3 n at the particle singularity. In fact by using the behavior of the solutions 
eq.([L3|) at the particle singularities, one obtains 

P nz = ^^ + P nz g'(z n ). (22) 

The fact that the time evolution is governed by the function g(z) will be the common 
denominator of the following developments. 



7 Conservation laws 

For a fuchsian differential equation of type (|13|) the following Fuchs relations hold p3| 



= £(l-/4 + 2/3 n <)+£(-3 + 2/3 B z B ). (23) 

n B 

The former expresses the fact that the function Q(z) behaves as 1/z 2 at infinity (Fuchs 
condition); the latter provides the second order residue of Q(z) at infinity (the total mass of 
the system) in terms of the other parameters. 

Below, we shall shows that these conditions implie some non-trivial conservation laws. 
We start recalling the equation of motion for tt z z 

7r 5 2 = 2e 2a d z (e- 2cT d z N) + 2n\d z N z + N*d z i? z (24) 

which by means of the eqs. fll~8|) and ( p!9[) can be recasted as follows 

n\ = -4d M Q(z) + d z 7i 2 z g + 2v\d z g. (25) 
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We have used a standard result about the Liouville equation that permits to connect a to 

Q(z) 

Q(z) = ~[dl(2a) + i(0,(2*)) a ]. (26) 

In a field theoretical language Q(z) is proportional to the energy momentum tensor associated 
to the Liouville equation. 

If we integrate both sides of eq. (53) over a closed path in the complex plane that encircles 
all the singularities, we find 

-§E P - = ^d M (£ Pn + E Pb) ~ (27) 

dt n n B 



2ti f dz (d z ii z z g + 2n z z d z 



Since both the sum over momenta and that over the first order residues vanish, the remaining 
integral must vanish as well. Recalling that tt z z ~ 1 / z 2 for large radius, we must conclude 
that g can grow at most linearly at infinity. Moreover, if we combine this information with 
the fact that the only singularities of g at finite radius can be poles, we find that g must be 
a rational function. 

Next, we multiply eq. (|25|) by z and we integrate again over a closed path in the complex 
plane to obtain 

~77 ^ ] Z'nPnz (28) 

at n 

47T<9 M (E Z nPn + E Z bPb) = 1 _ J~- 
n B ^ 

The previous equation combined with its complex conjugate gives both the law of conserva- 
tion of angular momentum 

L = Ti E K*nPn, ~ ZnPnz) = (29) 

at n 

and the "generalized conservation law" 

J t Y.( Z nPnz + Z n P nS ) = 2(1 - ^). (30) 
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The fact that the terms depending on g in eq. (|25| ) do not contribute can be easily seen if we 
write the r.h.s of (^) times z in the form 

- 4zd M Q(z) + 7i\{zd z g -g) + d z (z7r%g). (31) 

The last term is harmless since it is a total derivative of a rational function; the previous 
one falls again as 1/z 2 since the combination (zd z g — g) cancels the linear behavior of g at 
infinity. 

Eq.([ZBD is a severe constraint on the form of the reduced hamiltonian 7i describing the 
motion of the particles when gravity is eliminated. In fact J2 n Pnz = imposes 

H = H(( 2 ,...(n,P2,..-Pm) (32) 
with ( n = z n — Z\ and eq. (|28|) imposes that 

The general solution of this equation has the form 

H = £ (log(P„ 2 P n2 ) + log(C„Cn) a ' 1 ) + n (34) 

rt>l 

where 7i is a homogeneous function of degree zero in the variables P nz and l/£ n and Z)n>i(l — 
a n ) — 1 — ¥-. This property will be seen explicitly realized for two particles with TLq = 0. 

8 The two body problem 

In the two body problem, ir z z has no zero i.e. we are in absence of apparent singularities; 
g(z) = z/P\ z (z2 — z\) in order to have a non rotating frame at infinity. The simple residues 
/?2 are provided by the Fuchs relations in terms of particle masses mi,m 2 and the total 
energy. / appearing in eq.flTBD is given by the ratio of two hypergeometric functions which 
provides the metric || [7|, IHfl . On the other hand it is interesting that the particle trajectories 
can be determined even before the explicit computation of the metric. In fact we have for 
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the two particles the equation of motion ( |20| , |2~T|) . The rewriting of the two formulas for the 
problem at hand gives for I — z\ — and P z = P% z = —P<i z 

t-±. A— f (35) 

whose general solution is 

1 = const [(1 -fi){t- t ) -iL/2]^. (36) 

The scattering angle is immediately computed from the previous equation to be 9 scat t = 
irM/ (4tt — M) |l], EJ] . It is interesting to note that the scattering angle depends only on the 
total energy and not on the masses of the individual particles. The two equations (|35|) are 
generated by the hamiltonian H e ff = ln(Pj + Py) + /xln(^ + ly) which is unique up to an 
additive constant. Such an hamiltonian should not be confused with the total energy of the 
system; it is the effective hamiltonian for the time development of the two particle system 
and contains as an essential parameter the total energy M = 47r/x. 

9 The Af body problem 

We shall consider in this section the A/"-particle case with particular attention to Af = 3. 
In flat space the number of Lorentz invariants is easily seen to be 3A/" — 3 as for n > 3 
the 2 + 1-vector P n is determined by the scalar product of P n with Pi,P2,P3- In curved 
space we have the same number of invariants. As it is well known here the invariants are 
replaced by the traces of the holonomies around the world lines of an arbitrary collection of 
particles Jl], |25|, EJJ. The parallel transport monodromy matrices which belong to 50(2, 1) 



or in the fundamental representation to SU (1, 1) are defined up to a conjugation. Keeping in 
mind that every 577(1, 1) holonomy has three real degrees of freedom we have ?>N degrees of 
freedom to which we have to subtract the three degrees of freedom of the 577(1, 1) conjugation 
thus reaching as expected the same number of invariants as in the flat case. 

In the present approach parameterizing the dynamics in terms of the position of the par- 
ticles and their momenta is not a convenient choice, because these, as we have anticipated, 
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are not the natural variables for the associated Riemann-Hilbert problem. In order to un- 
derstand the problem better let us consider the case M = 3. It is useful to go over from the 
equation in projectively canonical form (|TB|) to the equivalent one presented below [E3 



y" + py' + qy = o 



with 



Q 



p p 
q ~2 ~T 



(37) 



(3c' 



The transformation can be so chosen as to have the following Riemann scheme 

/ 1 z 3 z A oo i 
Poo 

\ Ml M2 A*3 2 Poo + /loo J 

where the Fuchs relation 

Ml + P2 + A*3 + Moo + 2poo = 1 

fixes Poo in terms of the other parameters and 



(39) 



(40) 



p(z) 



M2 1 - M3 



1 



2 — 1 -2 — -23 -2 — -2A 
* 3 (*3-l)# 



z(z-l) z(z - l)(z - z 3 ) 
z A (z A - l)b A 



z(z - l)(z - Z A ) 

with k = Poo(Poo + Poo)- Here we use the unconstrained parameterization of ref.j23|; the 
connection with the parameters p„ and (3 n can be easily obtained by comparing the residues 



in eq.(|38 



The parameters p^ and p^ are free real parameters which give the trace of the holonomies 
around the particles and infinity (total energy). z A is given in terms of the z n and the P n 
(n = 1,2,3), while H is given by the no-logarithm condition [] Q A (z A ) = (3 A /2, i.e. the 



1 Here Qa(z) denotes the function Q without its poles in z A - 
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requirement that the point za is a regular point for 2a Such H is easily found 



H = Z A {Z A - 1){Z A - Z 3 ) | 6 2 _ + 

23(^3 - 1) I A 

+ + (41 ) 

- 1 - ^3 2U(^4 - 1) J 

Thus we have the following free real parameters ^1,^2, A*3, /Axu Re(^A), Ihi(za)j Re(£>A), im (&A) 
which are eight real parameters, to be contrasted with the six real invariants we have enumer- 
ated at the beginning of this section. The reason is that the constraint that the holonomies 
described by eq. fl3?D belong to 577(1,1), fixes the value of Re(&A); Ihi(&a) i n terms of all 
other parameters. The explicit form of this relation has been investigated in the mathemat- 
ical literature J28| but we are not aware of any explicit form of it. Adding a new particle 
introduces a new particle singularity with exponent (14 i.e. p(z) acquires a new term of the 
type (1 — fJL^)/(z — £4), but in addition a new apparent singularity is generated say at z B . 
The indices of this new singularity have to be integers, and 2 while the new simple residue 
&b has to be a function of all other parameters, and in particular of Zb, because the SU(1, 1) 
character of the holonomies allow only an increase of 3 in the number of free real parameters 
i.e. /i 4 , Rezs, Imzs- Adding other particles does not alter the procedure. 

The complete solution of the dynamical problem requires that we provide the equations 
of motion for all the variables parameterizing the dynamics of the system and in particular 
for the apparent singularities and their first order residues. To this purpose, we shall follow 
here a synthetic approach which will encode in a single equation all the information on the 
time development of the positions and residues of all singularities. To accomplish this we 
shall exploit the trace part of equation for i.e. 

d(2a) 



dt 



N z d z {2a) + d z N z + N z d- Z {2a) + d- z N z (42) 



which combined with the time development of tt z z (z), eq.([24]), provides the time derivative 
of the reduced conformal factor a i.e. 

2a = g{z)d z {2a)-g'{z) + C .c. (43) 
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We shall exploit now the fundamental relation ( p6|) of the theory of fuchsian differential 
equation [[23|, |2~T| . Taking the time derivative of eq.(^) we have 

Q(z) = \g"\z) + 2g'(z)Q(z) + g{z)Q'{z). (44) 

The above equation implies that the time development of Q(z) is given by a conformal 
diffeomorphism generated by Q(z) itself. We are not in presence of a classical transformation 
but of an anomalous one. Extensions of the classical Liouville theory have been considered 
27J. 



in 



It is interesting that this equation, in addition to the particle equations of motion, con- 
tains the whole Gamier system for the time evolution of the apparent singularities. We recall 
that g(z) is the rational function introduced in eq. (f21~D ; it is given by a sum of simple poles 
located at the apparent singularities and of a first order polynomial. 

Equating the simple residue of eq. fl44]) at infinity one simply obtains 



|(E^ + E/5b)=0 (45) 

which is already contained in the Fuchs relation (^), while equating the second order residues 
at infinity we obtain 



d 
If 







n n B 

which using the second Fuchs relation (^3|) provides d^/dt = 0, where = 1 — M/An i.e. 
the conservation of total energy. On the other hand if we equate the second order residues of 
eq.(|4]) at the particle pole z n we obtain d\j?Jdt = 0, where \i n = m n /47r i.e. the conservation 
of the individual particle rest masses. Matching of the third order residue on the particle 
singularity z n gives 

Z n = -g{z n ) (46) 

while matching of the first order residue gives f3 n i.e. 

4 = ^V(^n) + g\z n )L3 n . (47) 
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The motion of the apparent singularities is given by the matching of the residues of the third 
order pole at za- If we Laurent expand g(z) around za 

g(z) = + g + gi (z- z A ) + 0((z - z A f) (48) 

z- z A 

we obtain 

za = -g-ipA - 9o (49) 

while the first order pole gives 

Pa = -2g_ x Q' A (z A ) + giP A ~ \q2- (50) 

Equations (|49|) and ( |50|) can be written in hamiltonian form; for example in the three body 
problem, by choosing g(z) as to cancel the pole in eq.(|TP|) and keep particle 1 fixed at Z\ — 
and particle 2 fixed at z<i = 1, we have 

dz A z A {z A - l)(z A - z 3 ) 



dz 3 z 3 (z 3 - 1 

p.- 1 1 



za 1 z A 



dH , . 



and 

^ = (52) 
dz 3 dz A ' 

where H is the so called Gamier hamiltonian and it given by (f4~T|). The parameter (3a is 
connected to those appearing in H by the relation with 

£ = + + + (53) 

2 2\ z A za - 1 z a z 3 J 

The explicit expression of H can also be written for the general M body problem. 

It is now possible to interpret the problem of gravity in presence of particles a la Einstein- 

Infeld- Hoffmann [E9| i.e. as pure gravity without particles in which one allows for singular 

solutions of the field equations. What we found is that the deficit of the singularity has to 

remain constant (equivalent to the constancy of rest masses) , the total angular deficit has to 

remain constant (equivalent to the constancy of the total energy) and the field singularities 

follow geodesies of space- time. 
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10 Conclusions 



The application of the instantaneous York gauge K — 0, to 2 + 1 dimensional gravity 
is practically restricted to the case of open universes; on the other hand in this case the 
advantages are several and profound. The momenta n a b conjugate to the space metric have 
only two independent components and the diffeomorphism constraint has a very simple 
solution i.e. the n z z (n z 2 ) are rational functions of z (z) whose residues are the particle 
momenta P n . 

The advantage is still greater with regard to the hamiltonian constraint which in this 
gauge is equivalent to an inhomogeneous Liouville equation for the reduced conformal factor 
a; the sources are the particle singularities i.e. the poles of n z z and the apparent singularities, 
i.e. the zeros of n z z . The function N is given by the derivative of a with respect to the total 
energy while N z is given in terms of N by 

N z = ~d x N + g(z), (54) 

7T 

where g(z) is a rational function which cancels the polar singularities of the first term and 
grows at infinity not faster than z. The linear term in g(z) is fixed if we want to deal with 
a reference frame which does not rotate at infinity. 

The equations of motion for the position of the particles and of the apparent singularities 
and also the time dependence of the linear residues at such singularities are given by the time 
development of the function Q(z) entering the fuchsian differential equation which underlies 
the conformal factor describing the space metric. Such a time development is given by the 
energy momentum tensor of a conformal Liouville theory and it encodes completely the 
dynamics of the system. 
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